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We investigate the discovery potential of the LHC experiments for i?-parity violating supersym- 
metric models with a stau as the lightest supersymmetric particle (LSP) in the framework of minimal 
supergravity. We classify the final states according to their phenomenology for different i?-parity 
violating decays of the LSP. We then develop event selection cuts for a specific benchmark scenario 
with promising signatures for the first beyond the Standard Model discoveries at the LHC. For the 
first time in this model, we perform a detailed signal over background analysis. We use fast detector 
simulations to estimate the discovery significance taking the most important Standard Model back- 
grounds into account. Assuming an integrated luminosity of 1 fb _1 at a center-of-mass energy of 
^J~s = 7 TeV, we perform scans in the parameter space around the benchmark scenario we consider. 
We then study the feasibility to estimate the mass of the stau-LSP. We briefly discuss difficulties, 
which arise in the identification of hadronic tau decays due to small tau momenta and large particle 
multiplicities in our scenarios. 
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INTRODUCTION 



One of the main objectives of the experiments at 
the Large Hadron Collider (LHC) is to search for new 
phenomena beyond the Standard Model of particle 
physics (BSM) at and above the TeV energy scale. A 
well motivated scenario is the i?-parity violating min- 
imal supergravity model (mSUGRA) with baryon tri- 
ality (B3) |l|-|5|. Contrary to the i?-parity conserving 
mSUGRA model (also called the constrained MSSM) 
IE 0i lepton number is violated. This leads to natu- 
rally light neutrino masses without introducing either 
a new see-saw energy scale at My = O(10 10 GeV) 
or new gauge-sing let superfields [| i-[nj]. However, 
within this model the lightest supersymmetric parti- 
cle (LSP) is not stable. Thus the lightest neutralino is 
no longer a viable dark matter particle [9(J and other 
candidates must be considered, such as the axino 
16l . gravitino fr^Jlfjj , the lightest U -parity particle 
HO, III or slightly modified models like the NMSSM 
with a gravitino LSP (22J. More importantly for the 
focus of this paper, in this case, the LSP is no longer 
bounded by cosmology to be the lightest neutralino 
23 1 . Any other supersymmetric particle is possible 
24 , [25| . This can lead to dramatically different sig- 
natures at colliders HiS. Within the B 3 mSUGRA 
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model, certain non-neutralino LSP candidates are pre- 
ferred [2^, [29|, H3] ■ In particular for small i?-parity vi- 
olating couplings the lightest scalar tau, the stau, is a 
possible LSP in large regions of the mSUGRA param- 
eter space. This has been known for the i?-parity con- 
serving case for a long time, but has been discarded for 
the above cosmological reasons. In the i?-parity violat- 
ing case the stau is a natural LSP. We note however 
that the stau can also be the lightest supersymmet- 
ric particle within the (i?-parity conserving) MSSM 
spectrum if one adds another (lighter) particle to the 
spectrum, like the gravitino. The stau then decays 
into this new particle and can be again consistent with 
cosmological observations; see e.g. Ref. (3lj . 

In this paper, we consider in detail the discovery 
potential at the LHC for the B3 mSUGRA model with 
a stau LSP. We focus on the case of early data at 
7 TeV center-of-mass energy. In order to be specific, 
we mainly restrict ourselves to one of the B3 mSUGRA 
benchmark points discussed in Ref. [32|], as well as 
related scenarios. As we shall see, besides the good 
theoretical motivation, the stau LSP scenarios are also 
readily testable with early LHC data, because they can 
lead to a multi-lepton signature that is hard to achieve 
with the SM interactions. 

The outline of the paper is as follows. In Sec. HH wc 
present the model, which we investigate and then dis- 
cuss in detail how a stau LSP can arise. We also review 
the benchmark scenario that is relevant for this work. 
An overview of all expected stau LSP signatures at the 
LHC is given in Sec. IIII) i.e. we consider the differ- 
ent dominating i?-parity violating operators. We then 
present in Sec. IIVI in detail the particle multiplicities 
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and kinematic properties of our benchmark scenario 
and for the most important SM backgrounds. Based 
on this, we develop in Sec.|V]a set of cuts that allow a 
discovery of these stau LSP scenarios with early LHC 
data. We employ a fast detector simulation and esti- 
mate systematic uncertainties in our analysis. We also 
investigate the possibility of reconstructing the stau 
LSP mass. Finally, we comment in Sec. IVII on poten- 
tial difficulties with tau identification due to the large 
particle multiplicities in our scenarios. We summarize 
and conclude in Sec. IVII1 

In App. [XJ we review some basic properties of the 
investigated benchmark scenario. We give an overview 
over different definitions of significances in App. [Bj 



II. THE 7?-PARITY VIOLATING MSUGRA 
MODEL 

A. Motivation 

In the B 3 mSUGRA model HHH] there are six free 
parameters at the unification scale (Mqut) 



M , Mi/2, A), tan/?, sgn(/t), A. 



(1) 



Here Mo, Mi/ 2 are the universal supersymmetry 
breaking scalar and gaugino masses, respectively. A 
is the universal supersymmetry breaking scalar trilin- 
ear coupling and tan /3 is the ratio of the two vacuum 
expectation values. sgn(/i) = ±1 is the sign of the 
Higgs mixing parameter. The magnitude of y, as well 
as the respective bilinear scalar coupling Bq are fixed 
by radiative electroweak symmetry breaking [33j. In 
the B3 mSUGRA model the superpotential is extended 
beyond the MSSM by the following terms 



24] 



Wb 3 = XijkLiLjEk + X'ijkLiQjDk + KiLiH2 ■ (2) 

Here Li, Qi denote the lepton and quark SU(2) dou- 
blet superfields. H 2 denotes the Higgs SU(2) doublet 
superfield which couples to the up-like quarks. Ei, Di 
denote the lepton and quark SU(2) singlet superfields, 
respectively. i,j,k € {1,2,3} are generation indices. 



Xijk denote nine (anti-symmetric in i H j), X' i j k 
twenty-seven dimensionless couplings. Ki are three di- 
mensionful parameters, which vanish at the unification 
scale 0. All the operators in Eq. ([2]) violate lepton 
number. 

The parameter A in Eq. ([!} goes beyond the R- 
parity conserving mSUGRA model. It refers to a 
choice of exactly one of the thirty-six dimensionless 
couplings in Eq. ([2]) at a time. 

A e {Xi jk , X'i jk }, i,j,fc= 1,2,3. 



Given one coupling at the unification scale, through 
the renormalization group equations (RGEs) other 
couplings are generated at the weak-scale (Mew) [U 

Hi. 



In principle, one can also choose more than one im- 
parity violating coupling at the unification scale. We 
mainly restrict ourselves to only one coupling, because 
it makes the investigation of the parameter space more 
manageable. However, our approach is also motivated 
from a phenomenological point of view, because exper- 
imental constraints on products of different i?-parity 
violating couplings are in general more restrictive than 
on only a single coupling [12j. Furthermore, in the SM 
the top Yukawa is at least an order of magnitude larger 
than the other Yukawa couplings. 



B. Stau LSP 

As mentioned above, in the B3 mSUGRA model the 
LSP is no longer constrained to be the lightest neu- 
tralino, X\. Since the RGEs are modified compared 
to the i?-parity conserving mSUGRA model by the 
interactions in Eq. @, as well as the corresponding 
soft supersymmetry breaking terms, in principle new 
LSPs may arise as a function of the input parameters 
in Eq. ([!]). For large values of specific couplings, i.e. 



A ~ O(10 2 ) , this does in fact happen 

p, R for large Ai 32 , 

Vi for large X' ijk , 

e R for large Ai2i,Ai3i,A 23 i • 



(4) 



as discussed in detail in Ref. [26|, |29|, |30|. The rele- 
vant couplings are given here in parentheses. How- 
ever, even for small values of arbitrary couplings there 
are large regions of parameter space, where the light- 
est stau, fi, is the LSP. This occurs for large values 
of Mi/2 (which drives up the X? mass faster than the 
stau mass) and small values of Mo; see App. I A 21 for 
the tan j3 dependence. For a fixed value of Mo there 
is always a value of M x / 2 above which the stau is the 
LSP For example, at A = -100 GcV, tan /? = 10 
and sgn(/i) = +1 the region is approximately given by 



LSP 



M 



1/2 



3.8 • M + 175 GeV. 



(5) 



It is the purpose of this paper to study in detail the 
detectability of stau LSP models with early LHC data. 

There is a series of papers in the literature on a 
stau LSP @, M, S3, HI [33-11} • To the best of our 
knowledge a stau LSP was first considered in Ref. (37$ 
in a bilinear imparity breaking model with a focus on 
charged Higgs phenomenology [9l|. There has been 
further work on a stau LSP in bilinear imparity violat- 
ing models [II SI SI- In Ref- 0, LHC phenomenol- 
ogy is considered in some detail, however the focus is 
(3) on a Xi LSP, with only peripheral mention of the stau 



LSP discovery reach. In Ref. [40j the decay lengths 
of slepton LSPs arc considered in both bilinear and 
trilincar i?-parity violating models. LEP II searches 
for a stau LSP were performed in Refs [45l - |47j . allow- 
ing only for two-body imparity violating decays of the 
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BC1 


BC2 


BC3 


Operator 




X311L3Q1D1 


A331L3Q3 Di 


Coupling 

{@M G ut) 


0.032 


3.5 • 10" 7 


0.122 


Coupling 

(@M EW ) 


0.048 


1.1 ■ 10~ 6 


0.344 


LSP 


fi 


n 


v T 


M (GeV) 








100 


Mi/a (GeV) 


400 


400 


250 


tan P 


13 


13 


10 


A (GeV) 








-100 



TABLE I: The main properties of the B 3 mSUGRA bench- 
mark points. 



stau. Recently, Ref. [43| investigated stau LSP sce- 
narios, where the stau decays leptonically via 2-body 
decays. It was shown that these scenarios might al- 
ready be discovered in early LHC data via same-sign 
trilcpton events. 

In Ref. Q the B3 mSUGRA model was constructed. 
It was shown that extensive regions of parameter space 
lead to a stau LSP. Furthermore it was shown that in 
some regions of parameter space the stau dominantly 
decays via 2-body decays and in others via 4-body 
decays. In Ref. [32J benchmarks were defined for ex- 
emplary phenomenological and experimental analyses, 
for both the 2-body and 4-body decays. These are 
discussed below. In Ref. [27} , the origin of the 2-body 
and 4-body decays from the RGEs was studied in de- 
tail. The example of resonant slepton production was 
then analyzed in the case of a stau LSP. 

We here present for the first time a comprehensive 
analysis at the LHC of stau LSP scenarios. We focus 
mainly on one specific benchmark point BC1 [32| . For 
this we investigate the signal and the background in 
great detail. We include a fast detector simulation 
using the program package Delphes 50]. 



C. Benchmark Scenarios 

In Ref. [HI four Bonn-Cambridge (BC) benchmark 
points for studying i?-parity violating models in de- 
tail were defined. BC1-BC3 are B 3 mSUGRA mod- 
els, BC4 involves a baryon number violating opera- 
tor. The main properties of the points BC1-BC3 are 
summarized in Table fl] BC1 and BC2 have a stau 
LSP, the focus of our investigation here. In contrast 
to BC2, BC1 involves an additional purely leptonic 



operator, \\1\L\L1E\. Thus the Xi and stau decay 
purely leptonically. In the production and cascade de- 
cay of squarks and gluinos at the LHC this leads to 
many leptons in the final state. In BC2 the stau LSP 
decays purely hadronically, leading to significantly less 
final state leptons. Here we focus on the BC1 scenario, 
due to the better prospects of early discovery at the 
LHC. The mass spectrum as well as the dominant de- 
cay branching fractions for BC1 are given in App. [Al 



III. STAU LSP SIGNATURES AT THE LHC 

As a first step of our analysis, we classify the 
main LHC signatures, assuming that the stau decay 
is dominated by only one i?-parity violating operator, 
cf. Eq. ([2]). For simplicity we focus on the cascade 
process 



11/99 -> 11 



JJTTT1T1 



(6) 



where q is a squark, and j denotes a jet. The final- 
state staus can only decay via i?-parity violating op- 
erators. We then classify the signatures according to 
the possible stau decays. The results arc summarized 
in Tables HI HH and 



LiQjD k , or UiDiD k 



IVj for the operators LiLjEk, 
92|, respectively. Here we as- 
sume that always only one operator is dominant (93[. 
In the left column we denote the dominant coupling, 
in the middle column, the corresponding dominant 2- 
or 4-body stau decays. The resulting final state signa- 
tures for the cascade of Eq. are given in the right 
column. Here £ denotes an electron or muon, r a tau 
and ~$t missing transverse energy due to neutrinos in 
the final state. We have not included neutrinos from 
tau decays. 

We show in Fig. [T] as an example the 4-body stau 
LSP decay into t + fi + e~D e via Am, cf. the first decay 
in Tab. |TTJ The stau does not directly couple to the 
LiL 2 Ei operator and thus first couples to a virtual 
neutralino radiating off a tau lepton. The neutralino 
then couples to, for example, an electron anti-neutrino 
and a virtual electron sneutrino, 9 e , which decays via 
A121 into an electron and a muon. We end up with a 
4-body decay of the stau LSP. 

We observe that the partial width, T4, for the 4- 
body decay gets suppressed if the sfermion (gaugino) 
masses increase, i.e. the behavior is approximately 
T4 ~ m^ 4 (m~o) [||. Furthermore, we see that for 
the scenario BC1, the decay is mainly mediated by 
the Xi- The (mainly right-handed stau LSP) couples 
much stronger to the (bino-like) Xi than to the heavier 
(wino- and higgsino-like) neutralinos. If we also take 
into account the Majorana nature of the neutralinos 
we obtain from Fig.[T]the second stau LSP decay mode 
of Tab. HU 

The stau LSP can in principle also decay via a vir- 
tual chargino instead of a virtual neutralino. However, 
these decays are suppressed by several orders of mag- 
nitude. On the one hand, they are suppressed due to 
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FIG. 1: Feynman diagrams contributing to the 4-body stau LSP decay fj — > t + n + e~ v e via Am. In this example, the 
stau LSP decays via a virtual neutralino Xi (J = 1 . . .4). In addition, the first, second and third diagram involve a virtual 
electron sneutrino, v e , a virtual (left-handed) smuon, fiL, and a virtual right-handed selectron, Sr, respectively. 



coupling 


fx decay 


LHC signature 


Al21 = — A211 


r + n + e~v e 
r + fi~e + u e 
r + e + e~P M 
T + e~e + v^ 


2j + At + U + $ T 
with I = e, [i 


Al22 = — A212 




A131 = — A311 


e + v e 


2j + 2r + 21 + $ T 


A132 = — A312 




A231 = — A32I 


e + L' li 


A232 = — A322 




Al23 = — A213 




A133 = — A313 


e + i> T 

T + V e 
T V e 


2j + 2t + 21 + $ T 
2j + 3r + U + $ T 

2j + 4t + $ t 


A233 = — A323 





TABLE II: Scenarios assuming one non- vanishing LiLjE k 
operator. The column in the middle shows the possible 
stau LSP decays, for the Xij k couplings in the left column. 
The right column shows the resulting LHC signatures. The 
SUSY cascade qq/gg — > — >■ 33X1X1 3JTTTif\ has been 
assumed. Note that gluino, g, pair production instead of 
squark pair production will usually give two additional jets; 
for example via the decay g — >• jq. 



heavy propagators. On the other hand, the (mainly 
right-handed) stau couples in mSUGRA more strongly 
to the (bino-like) Xi than to the (wino-like) xt ■ 

Note, that due to the Majorana nature of the Xi, 
every charge combination of the two staus is possible. 
The flavor and charge of the leptons and quarks in 
the final state, are determined via the different decay 
modes of the stau LSP. 

In general, more complicated decay chains than 
Eq. ([6]) can occur, leading typically to additional fi- 



coupling 


fx decay 


LHC signature 


Kjk 


T + iijdke + 
r + Ujdke~ 
T + djd k v e 
T + djd k v e 


6j + 4t + ££ 
6 3 + 4r + i + $ T 

63 + 4r + f, T 




T + Ujdk[l + 
T + Ujdk/J,~ 

T + djdkVfj. 
r + djd k u^ 


A 3j fc 


Ujd k 


6j + 2t 


TABLE III: Same as Tab. [TT] but for one non-vanishing 
LiQjDk operator. 


coupling 


f + decay 


LHC signature 




T + mdjdk 
r + Uidjd k 


8j + 2r 



TABLE IV: Same as Tab. HQ but for one non-vanishing 
UiDjD k operator. 



rial state particles. The most important are: 

• Additional jets from the production of gluinos 
and their subsequent decays into squarks and 
quarks. 

• Additional leptons from decays involving heavy 
neutralinos and charginos. For example, a left- 
handed squark might decay to a X2 or xt > which 
then decays to a lepton and a slepton. 



• Additional leptons from the decay of a Xi next 
to-next-to-next-to LSP (NNNLSP) into right 
handed sclcctrons or smuons 



i.e. Xi — > 



followed by 



The last process is special for stau LSP scenarios. 
Within mSUGRA, it is kinematically only allowed if 
M 1/2 > M , e.g. for M 1/2 > 400 GcV for M = 
GeV. Note, that M1/2 increases the mass of the (bino- 



like) Xi faster than the mass of the £r [5l|. If these 
new decay channels are open, we can have two addi- 
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tional taus and four additional electrons or muons in 
the final state. See Ref. [13] for explicit examples. 

The mufti charged lepton final states (especially 
electrons and muons) are the most promising signa- 
tures to be tested with early LHC data. It is relatively 
easy to identify electrons and muons in the detector 
and for high multiplicities the SM background is very 
low [56|. In the next section, we therefore investigate 
the discovery potential of stau LSP scenarios, where 
the stau decays in the 4-body decay mode via the op- 
erator \\2\L\L2E\, leading to a maximal number of 
muons and electrons in the final state, cf. Tab. HU We 
also analyze the potential of reconstructing the mass 
of the stau LSP, once a discovery has been made. 

Note, that scenarios, where the stau decays via a 2- 
body LiLjEk mode (i,j or k = 3) might have a similar 
discovery potential, because the charged leptons have 
on average larger momenta. We also expect larger 
values of fix, because a 2-body stau decay via LiLjEk 
involves always one neutrino, cf. Tab. [Til 

If the stau LSP couples mainly to a LiQjDk oper- 
ator, there are (at parton level) at least six jets and 
two taus in the final state, cf. Tab. IIIII Furthermore, 
for i = 1,2, there are possibly one or two additional 
charged electrons or muons as well as four additional 
taus. Again, due to the Majorana nature of the Xi, 
all charge combinations for these electrons and muons 
are possible, leading for example to like sign dilcp- 
ton events. Unlike the scenarios shown in Tab. |TTJ 
(LiLjEk), the scenarios in Tab. IIIII (LiQjDk) do not 
necessarily lead to fx (apart from neutrinos from tau 
decays). It should therefore be possible to directly re- 
construct the mass peaks of sparticles, especially the 
stau LSP. However, combinatorial backgrounds due to 
the many jets in the final state complicate this task. 

A special case arises if the LiQjDk operator in- 
volves quark superfields of the third generation, i.e. j 
or k = 3. For k = 3 two jets in the final states are 
6-jets, one from each stau decay. For i = 1,2 and 
j = 3 the situation is similar. We also obtain one b- 
jet from nearly each stau LSP decay. However, decays 
which involve an electron or muon are now kinemat- 
ically suppressed or forbidden due to a top quark in 
the final state, cf. Tab. IIIII Finally, for a non- vanishing 
coupling Ag 3fc , the stau might only decay via a 3-body 
mode and a virtual top-quark. See Ref. [231 f° r details. 

In this paper, we focus on stau LSP scenarios with 
a lepton-number violating LiLjEk operator. However, 
we briefly discuss signatures, which arise from a non- 
vanishing baryon- number violating operator UiDjDk- 

We see in Tab. [TV] that a non-vanishing \"j k cou- 
pling leads mainly to a 4-body decay of the stau LSP, 
resulting in at least eight parton level jets and four 
taus in the final state. Due to the many jets, we ex- 
pect tau identification via its hadronic decay modes to 
be very difficult; see Sec. IVII A jet from a cascade de- 
cay might overlap with the tau jet. Furthermore, the 



three jets in a stau decay can be boosted such that 
they might appear as only one jet. However, these 
jets might still be revealed by investigating the sub- 
structure of jets. See Ref. [52| for X? LSP scenarios, 
with purely hadronic (UiDjDk) LSP decays. See also 
Rcfs. [5^M37| for related work. We expect that similar 
techniques will work for stau LSP scenarios. 

For the special case of j or k = 3, the baryon- 
number violating stau decays lead to two fe-jcts in the 
final state. For i = 3, the final state up-type quark in 
Tab. IIVI will be a top-quark. If the decay into a top 
quark is kinematically forbidden, the stau LSP might 
decay in a 5-body decay via a virtual top. 



IV. SIMULATION OF SIGNAL AND 
BACKGROUND 

In this section we perform a full Monte Carlo anal- 
ysis of the stau LSP benchmark scenario BC1 (with 
A121 = 0.032 at Mqut', cf. Tab. [TJ and the dominant 
SM backgrounds. Our signal process is pair produc- 
tion of all SUSY particles. We also employ fast de- 
tector simulations. The mass spectrum and branching 
ratios (BRs) of BC1 are given in App. [S] 



A. Major Backgrounds 

In what follows, we only consider SM backgrounds 
that can lead to at least one (parton level) electron, 
muon or tau in the final state. Furthermore, we ex- 
pect from most of the SUSY (signal) events additional 
energy in the form of hard jets, that arise from decays 
in the upper parts of the decay chain. 

We thus consider the following SM processes as the 
major backgrounds in our analysis: 

• ti production. 

• Z + jets. The Z can decay into a pair of charged 
leptons. Because the SUSY events posses addi- 
tional large amounts of energy from jets, we only 
consider Z production with at least one hard jet 
at parton level. 

• IU+jets. The W can decay into a charged lepton 
and a neutrino. In order for the W production to 
be competitive with the SUSY processes, we de- 
mand at least two additional hard jets at parton 
level 

• Di-boson (WW, WZ and ZZ) production [Hj]. 

Table [V] gives an overview of the background sam- 
ples used in our analysis. QCD di- and multi-jet events 
have been neglected. It has been shown e.g. in Ref. foOj 
that QCD background can efficiently be suppressed in 
multi- lepton final states. 
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Sample 


sub-sample 


simulated events 


Generator 


comments 


tt 




240 000 


MC@NL0 + Jimmy/Herwig 




Z + jets 


— — x z — : 

Z -> e + e + > 1 jet 
Z — > (i + > 1 jet 
Z ->• tt+ > 1 jet 


222 000 

231 000 

232 000 


Alpgen + Jimmy/Herwig 


for each process split in the 

excfusive sampies f , . . . , 4 extra partons 

and the inciusive sample > 5 extra partons 


W + jets 


W -> ez/ + > 2 jets 
W ->■ /w/+ > 2 jets 
W ->• ri/+ > 2 jets 


518 000 
642 000 
659 000 


Alpgen + Jimmy/Herwig 


for each process split in the 

exclusive samples 2, . . . , 4 extra partons 

and the inclusive sample > 5 extra partons 


di-boson 


WW 

WZ 

ZZ 


30 000 
20 000 
9 990 


Jimmy/Herwig 




signal 




«10 000 each 


Jimmy/Herwig 





TABLE V: Monte Carlo samples of the SM background and signal events used for our analysis. The third column shows 
the number of simulated events (for the Alpgen samples after MLM matching [7ll - [73| ) and the fourth column shows the 
employed generators. For the simulation of the signal we employed a special version of Herwig which also incorporates 
the 4-body decays of the stau LSP [g^]. Note, that we simulated approximately 10 000 signal events for each mSUGRA 
parameter point in Sec. IV Bl including BC1. 



B. Simulation and Selection Cuts 

The SUSY mass spectra were calculated with 
S0FTSUSY3 . [H|62[. The S0FTSUSY output was fed 
into ISAWIG1 . 200 and ISAJET7 . 75 [Hi in order to cal- 
culate the decay widths of the SUSY particles (beside 
the stau LSP). The signal processes, i.e. pair produc- 
tion of SUSY particles, was simulated with a modified 
version of Herwig6.510 [64 66] which also simulates 
the 4-body decays of the stau LSP [67j ■ We employed 
Jimmy4.31 [68[ to simulate the underlying event. 

The tt background was simulated with MC@NL03 . 41 
[6!l[7(j , the di-boson background in Herwig6 .510. We 
used Alpgen2 . 13 f7ll - l73j interfaced with Herwig6 . 510 
to simulate the Z+jets and TU+jets backgrounds. We 
also employed Jimmy4.31 in all cases for modeling of 
the underlying event. All Monte Carlo samples are 
based on the parton distribution functions given by 
CTEQ6L1 [74|. The employed Monte Carlo generators 



the ATLAS collaboration [60| for the expected back- 
ground uncertainties to calculate significances in the 
following. 



are also summarized in Table El 

Detector effects on signal and background were 
accounted for using the generic detector simulation 
Delphesl.8 Its detector settin gs w ere tuned to 
an ATLAS-like detector at the LHC [HI . The results 
of the detector simulation were cross checked with the 
PGS4 simulation [76| using its generic LHC tune. A 
sufficient agreement of the two codes was observed 
for most observables. However, the identification of 
electrons and especially tau leptons showed some dis- 
crepancies, cf. Sec. IVI1 Detector simulations with full 
detail of the calorimctry and tau identification algo- 
rithms need to be done to get more reliable results. In 
the comparison between PGS and Delphes we observe 
some shifts in the energy scale especially for electrons 
in our signal events. Using the differences between 
the two fast simulations as an estimate for the energy 
scale uncertainty would lead to unreasonably large sys- 
tematic uncertainties in the estimate of the expected 
discovery significance. Instead, we take an estimate of 



particle 


transverse momentum 


pseudorapidity 


electron 


p T > 7 GeV 


\n\ < 2.5 


muon 


P t > 6 GeV 


W <2.7 


tau 


Pt > 10 GeV 


M < 2.5 


jet 


p T > 20 GeV 


M < 5.0 



TABLE VI: Cuts for the particle selection for the signal 
and background. 

The particle selection was guided by the definitions 
used by the ATLAS collaboration for SUSY studies, 
cf. Ref. [6(| PP- 1518]. The selection cuts are given 
in Tab. IVII We similarly followed these guidelines for 
the overlap removal of reconstructed objects. This is 
needed, because a single particle may be reconstructed 
as several different objects. 

For the overlap removal we simply use the distance 
in pseudorapidity i] and azimuthal angle cf>, defined as 
AR = (A0) 2 + (A?y) 2 . Objects are selected in the 
following order: 

1. Muons, if no jet and no previously selected muon 
is present within AR < 0.4. 

2. Electrons, if no jet is present within 0.2 < AR < 
0.4 and no previously selected electron within 
AR < 0.4. 

3. Hadronically decayed taus, if no electron or tau 
has already been selected within AR < 0.4. 

4. Jets, if no electron, tau or another jet has already 
been selected within AR < 0.4. 

These cuts take care of the fact that electrons and taus 
are usually also reconstructed as jets. Furthermore, 
electrons are more reliably identifiable than taus. In 
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(d)Number of taus. 



FIG. 2: The number of reconstructed particles per event after selection cuts and overlap removal as described in Sec. lIV ill 
The color code for the curves in all four plots is given in (d). The different background contributions are stacked on 
top of each other, while the expected signal is shown in front of the histograms in red. The signal corresponds to the 
benchmark scenario BC1 [32l. The number of events are scaled to 1 fb" 1 at Js = 7 TeV. 



addition, one does not want to select electrons or 
muons, that stem from heavy flavor decays within jets. 



C. Particle Multiplicities and Kinematic 
Properties 

In this section we show the basic kinematic proper- 
ties of the BC1 scenario compared to the most impor- 
tant Standard Model backgrounds. We also motivate 
the cuts to obtain a good significance and a good signal 
over background ratio. No further event selection cuts 
beyond the object selection cuts given in Sec. IIVBI 
are applied here. All samples arc scaled to an inte- 
grated luminosity of 1 fb -1 at a center-of-mass energy 
of yfl = 7 TeV. 

We show in Fig. [2] the number of reconstructed jets, 
electrons, muons and taus, respectively The different 
background contributions are stacked on top of each 
other, while the expected BC1 signal is shown in front 
of the background histograms. Error bars correspond 



to statistical uncertainties of the generated samples. 
From the discussion in Sec. HIT] we expect 2-4 jets at 
parton level for the scena rio BC1. This behavior can 
also be seen in Fig. 2(a) (red histogram), where the 
signal distribution is maximal at 3 jets. In general, 
we can get in addition to the parton-level jets also 
jets from parton shower radiation and possibly from 
non-identified hadronically decaying taus. 



According to Fig. 2(a) the SM background is pre- 
dominant even for high jet multiplicities. The ti back- 
ground contributes via the two (parton-level) 6-jets 
from top decays as well as jets from a hadronically 
decaying W . Therefore, a 6-jet veto may in princi- 
ple suppress this background. However, the signal 
can contain 6-jets as well, and we will not consider 
6-tagging in our analysis. 



The bins in Fig. 2(a) with < 1 jet are dominated by 
the > 1 jet background (blue histogram). This is 
expected due t o the large Z + jets cross section. Note 
that Fig. 2(a) displays the event numbers on a loga- 
rithmic scale. In contrast, W + jets (green histogram) 
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(a,)pT distribution of the hardest jet. 
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(b)pT distribution of the 2nd hardest jet. 




50 100 150 200 250 300 350 400 450 500 

p W [GeV] 

(c)pt distribution of the 3rd hardest jet. 
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(d)px distribution of the 4th hardest jet. 



FIG. 3: pr distributions of the four hardest jets after object selection cuts and overlap removal. 



mostly dominates for > 2 jets, because we only show 
W production in association with two parton-levcl jets 



in Fig. 2(a) 



The transverse momentum, px, distributions of the 
four hardest jets are shown in Fig. [3] We observe 
in Fig. 3(a) that the pr distribution of the hardest 
signal jet is relatively flat over several hundreds of GeV 
whereas the distribution of the background falls off 
steeply. One can also see a peak in the pt distribution 
of the signal around 360 GeV. This peak is due to the 
mainly hard jets from the decay of a squark into the 
Xi. Note that the mass difference between most of the 
squarks and the X? in BCI is 400-500 GeV, cf. App.[H 
The invariant mass of the hardest jet with the decay 
products of the X? might thus allow a reconstruction 
of the squark masses. 



Gluinos decaying into a jet and a squark will in gen- 
eral produce softer jets, because the respective mass 
difference is only « 100 GeV. Therefore, these addi- 



tional jets are most probably found in Fig. 3(b) 3(d) 



We see in Fig. [3] that the signal over background 
ratio improves if one only allows for very hard jets. 
This is also reflected by the visible hadronic mass, 




1000 1200 1400 1600 1800 2000 

HT' [GeV] 

FIG. 4: The number of events per lfb -1 as a function 
of the scalar sum of the transverse momenta of the four 
hardest jets, HT' , after object selection cuts and overlap 
removal. 



HT' = J2j e t i-aPTi wu ich is the scalar sum of the 
transverse momenta of the four hardest jets. The re- 
spective distribution is given in Fig. 2] Note, that this 
cut alone (no actual cut is defined at this point) would 
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(a) px distribution of the hardest electron. (b)pr distribution of the 2nd hardest electron. 

FIG. 5: pr distributions of the two hardest electrons after object selection cuts and overlap removal. 
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(a)pT distribution of the hardest muon (b)pr distribution of the 2nd hardest muon 

FIG. 6: pr distributions of the two hardest muons after object selection cuts and overlap removal. 



still give an overwhelming amount of QCD background 
events. We thus need to also make use of the charged 
leptons in the final state. 

The most striking signature of the BC1 scenario is 
the large number of electrons and muons in the fi- 
nal state . Wc s how t he respective multiplicities in 
Fig. |2(b)] and Fig. 2(c) Even with an integrated lumi- 
nosity of only 1 fb~ x at y/s = 7 TeV, we expect several 
dozen events with four or more electrons in the final 
state! 

The signal distribution of the number of recon- 
structed electrons peaks at 3 whereas for the muons it 
peaks at 1. This is exactly the expected behavior from 
the parton-level signatures reviewed in Sec. IIII1 cf. es- 
pecially the row for Am in Tab. [TTJ The decay of two 
stau LSPs in a typical BC1 event leads at parton level 
to 2-4 electrons and 0-2 muons. Furthermore, addi- 
tional electrons or muons can arise from other particle 
decays; e.g. in BC1 from leptonic tau decays. 

It is interesting to note that the ratio of the num- 
ber of reconstructed electrons and muons carries infor- 



mation about the involved B3 coupling. Comparing 
Fig. |2(b)| with 2(c) we see that the ratio of the av- 
erage number of electrons to muons is roughly three. 
This is because Am couples two lepton superfields of 
the first generation to one of the second generation. 
Thus, stau LSP decays produce more electrons than 
muons, cf. Tab. [TTJ If we instead had Am 7^ 0, the 
situation would be reversed (modulo differences due 
to the reconstruction efficiencies). 



According to Fig. 2(b) it is possible to obtain a 
nearly background free sample by requiring more than 
four electrons in the final state! However, such a cut 
would also veto most of the signal events and is there- 
fore not well suited for a study of early data. Further- 
more, electrons can easily be faked by early photon 
conversions or low multiplicity jets. Therefore, we de- 
mand (as a cut in the next section) "only" the presence 
of at least one electron with pt > 32 GeV, another 
electron with pt > 7 GeV and at least one muon with 
Pt > 40 GeV in the final state. It is hard for the 
QCD background to fake an electron and a muon at 
the same time. In addition only the leptonically decay- 
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FIG. 7: The number of events per lfb" 1 as a function 
of the scalar sum of the transverse momenta of electrons 
and muons, ^2pt, after object selection cuts and overlap 
removal. 



ing Ws from tt production, leptonically decaying taus 
from Z — > tt production, and di-boson production can 
lead to an electron and a muon at parton level. We 
thus expect a strong suppression of the background 
from this cut, cf. Tab. ED 

We display the pt distributions of the two hardest 
electrons in Fig. [5] and those of the two hardest muons 
in Fig. |51 The pr distribution of the background falls 

we see 



off more rapidly than the signal. In Fig. 5(b) 



that for pt > 150 GcV the signal even dominates over 
the background. 

Due to their large multiplicity and large transverse 
momenta, leptons beyond the sub-leading electron and 
leading muon can contribute significantly to the en- 
ergy deposition of all leptons. The cut J2Pt > 
230 GeV, which we employ in the next section, ac- 
counts for the fact that the signal lepton ptS are on 
average larger than the background. The lepton mo- 
menta are large because they lie at the end of a decay 
chain of heavy SUSY particles, cf. Tab. IVTlIl Note, 
that J^Pt ls the px sum of all electrons and muons. 
The distribution is shown in Fig. [7] 

We show in Fig. |2(d)| the multiplicity of recon- 
structed (hadronically decaying) taus. The signal dis- 
tribution peaks at 0. At first glance this is surprising, 
because from Tab. [TT] we expect 4 taus in most of the 
BC1 SUSY events. 

However, with respect to tau identification (ID), the 
topologies in BC1 (cf. Sec. are special. Due to a 
large number of taus and jets from the SUSY decay 
chains, overlaps between different tau jets and other 
jets make standard tau ID via its hadronic decays 
very difficult. The low visible momentum of many 
tau leptons in this scenario complicates the tau lepton 
identification further. As we show in more detail in 
Sec. IVI1 tau ID in BC1 has an efficiency of at best 
20%-30%. The exact number depends on the working 
point on the cfficiency-vs-rejection curve of the tau ID 



algorithm, that is used. Furthermore, we observe a 
strong dependence on the nature of the fast detector 
simulation used, i.e. Delphes or PGS4 (cf. Fig. [17)) . In 
any case the ID efficiency is expected to be a factor of 
2-3 smaller in BC1, than e.g. in Z — > tt events, even 
for tau leptons of the same momentum. 

Although the number of parton-level taus is much 
larger in BC1 compared to the SM backgrounds, this 
docs no longer hold for the reconstructed taus. This 
can also be seen in Fig. [5J where we display the (visi- 
ble) pt of the two hardest identified taus. Even for 
large momenta, the background always exceeds the 
signal. Naively, we would expect the contrary. Like 
the electrons, Fig. [3J the taus result from the decay 
chains of heavy SUSY particles. In addition, the tau 
leptons in BC1 are mostly very soft, which reduces the 
ID efficiency further. Therefore, we do not employ the 
taus to improve the signal to background ratio. 

In Fig. [HI we present the missing transverse energy 
distribution for the signal and the backgrounds. Even 
though we investigate a B3 scenario, where the LSP 
does not escape detection, the missing transverse en- 
ergy, fix, can be significant in BC1. For example, 
a cut on the missing energy of fir > 400 GeV im- 
proves the signal over tt background ratio from ini- 
tially 0(0.001) to O(l), cf. Fig.lU 

The reason for this are the neutrinos from the stau 
LSP decays 



X? 



and the neutrinos from the successive r decays. How- 
ever, we will not explicitly cut on fix in the event 
selection in order to keep the analysis complementary 
to searches for i?-parity conserving SUSY. 

It should be kept in mind, that the discovery of a fir 
signal does not necessarily contradict ftp models, fix 
alone is not sufficient to distinguish i?-parity violating 
from i?-parity conserving SUSY. Further observables, 
like kinematic edges, will be needed to gain insights 
in the SUSY mass spectrum. Finally, a lepton (linear) 
collider may be needed to clarify the nature of the LSP 
and the source of the missing energy. 



V. DISCOVERY POTENTIAL AND MASS 
RECONSTRUCTION 

A. Cut Selection and Significances 

We now employ in addition to the pre-selection cuts 
on reconstructed objects (see Sec. lIVBl) further cuts in 
order to increase the signal significance and the signal 
over background ratio. Wc show, that with an inte- 
grated luminosity of only 200 pb" 1 at y/s = 7 TeV a 
discovery of the benchmark scenario BC1 is possible. 
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(a)pT distribution of the hardest tau 
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FIG. 8: (Visible) pr distributions of the two hardest identified taus after object selection cuts and overlap removal. 
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FIG. 9: Missing transverse energy, fir, distribution after 
object selection cuts and overlap removal. 



From the kinematic properties presented in 
Sec. IIV C[ we implement the following cuts: 

• px(lst fi^) > 40 GeV. We demand at least one 
muon with pr > 40 GeV in the final state. 

• px(lst e^) > 32 GeV. We demand an electron 
with pt > 32 GeV in the final state. 

• pr(2nd e*) > 7 GeV. We demand at least a 
second electron in each event. Note, that pt > 
7 GeV corresponds to the electron pre-selection 
cut, cf. Tab. ED 



• J2pt > 230 GeV. We demand the pr sum of all 
electrons and muons to be larger than 230 GeV. 

• HT > 200/300/400 GeV. We employ also dif- 
ferent cuts on the pt sum of the four hardest 
jets, namely 200 GeV, 300 GeV and 400 GeV, 
respectively. 

The cut selection was optimized with the help of 
the TMVA toolkit [73 . The cuts were chosen so as 



to improve the signal to background ratio, based on 
simulated annealing [78| . This was iterated with dif- 
ferent sets of event selection variables, while some cuts 
were left fixed. We found that slightly different cuts in- 
crease the signal significance marginally for our Monte 
Carlo samples, but only at the cost of a higher risk 
to accept QCD events, which we could not simulate. 
Several other variables, like the £ + £~ invariant mass 
as Z-veto or the W transverse mass, were tested as 
well. However, they turned out to be of less relevance 
in combination with the variables finally used. 



The cut flow is given in Tab. IVIII and visualized in 
Fig. [TUJ We also quote the signal to square root of 
background ratio, S/ \[B (second last column), and 
the significance Zq [79( (last column) assuming a rel- 
ative background uncertainty of 50%. For a definition 
of Zq and related measures of significance see App. [Bl 

We have chosen such a large systematic uncertainty 
for the SM backgrounds as a conservative approach, 
because uncertainties are expected from jet and lep- 
ton energy scales as well as the predictions of cross 
sections. The ATLAS collaboration for example esti- 
mates the systematic uncertainties for an integrated 
luminosity of lfb -1 to be 50% for the background 
from QCD multijet events and 20% for the background 
from tt, W+jets, Z+jets, and W pairs in the context 
of SUSY searches [60|. In our analysis we did not 
consider QCD multijet background explicitly. Addi- 
tionally the systematic uncertainties will be larger for 
the very first searches and we therefore use the more 
conservative estimate of 50% uncertainty. The uncer- 
tainties in Tab. IVIII stem from the limited statistics 
of the simulated events (see Tab. [Vj. Note that for 
VF+jets we only show events with at least two jets at 
parton-level 



Employing only selection 
Tab. EH 



cuts (second row in 
we observe that the number of SM back- 
ground events (sixth column) is several orders of mag- 
nitude larger than the number of signal events (sev- 
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FIG. 10: Cut flow in BC1 at ^/s = 7 TeV. As in 
Tab. IVIII entries are scaled to an integrated luminosity 
of jLdt = lftr 1 and error bars include statistical un- 
certainties only. The significance Zq, denned in App. [Bl 
assumes a background uncertainty of 50%. 



enth column) [96j]. Without any cuts, we expect at 
= 7 TeV roughly 300 signal events for an inte- 
grated luminosity of 1 fb . 

We obtain a great improvement of the signal to 
background ratio to 0(0. 1) by demanding one hard 
muon (pt > 40 GeV) and one hard electron (px > 
32 GeV) in the final state. At this stage the main 
backgrounds stem from leptonic decays of tt (second 
column) and from leptonic decays of di-bosons (fifth 
column). These processes can produce an electron 
and muon at parton level. In addition, we also get 
some contributions from the W^+jets and Z+jets back- 
grounds. Here, the second charged lepton of the first 
or second generation stems mainly from jets that are 
misidentificd charged lepton. After the first two 
cuts, we already obtain a S/V~B ratio of roughly three. 
However, if one takes into account systematic uncer- 
tainties, that are expected to be large at the early 
LHC run, one cannot observe a clear signal, cf. the 
last column in Tab. IVIII 

Demanding at least a second electron in the final 
state (fourth row) the situation further improves. We 
now have a signal to background ratio of roughly one 
and a S/y/B ratio of 8.4. Therefore, by demanding 
two (hard enough) electrons and one muon in the final 
state, a signal of the scenario BC1 might already be 
visible in early LHC data, even though the significance 
Zq is still small without further cuts. 

As we have shown in Figs. [5] and [SJ the electrons 
and muons of the signal events also have on average 
larger transverse momenta than the charged leptons 
from SM processes. We therefore demand in the sixth 
row in Tab. lVIll the pr sum of the electrons and muons 
to be larger than 230 GeV; cf. Fig. With this cut 
we obtain a number of signal events that is roughly six 
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times larger than the number of background events. 
A clear signal of BC1 should now be visible, because 
S/VB = 22 and Z = 4.9. 

Up to now, we have only employed the electrons and 
muons that stem mainly from the stau LSP decays. 
However, in a typical SUSY event at the LHC, we will 
first produce a pair of strongly interacting sparticles. 
These sparticles then cascade decay down to the LSP 
producing hard jets in the final state. As can be seen 
in Fig. |31 these jets are in general harder than the 
jets from SM backgrounds. Therefore, we demand in 
Tab. IVIII as the last cut the pr sum of the four hard- 
est jets, Fig. 21 to be larger than 200 GeV (third last 
column), 300 GeV (second last column), or 400 GeV 
(last column) , respectively. We show different cuts for 
HT' in order to show how the signal over background 
ratio improves with harder cuts on HT' . 

For HT' > 200 GeV, we already obtain a back- 
ground sample without di-boson events (fifth row). 
Because we only have produced a finite number of 
background events (see Tab. fV]) . we mark a back- 
ground free sample by "< 1.0". If we harden the cut 
further to HT' > 300 GeV, we also veto all Z+jets and 
W + jets events. Finally for HT' > 400 GeV, we ob- 
tain even a background free sample! Note that at the 
same time, the number of signal events is only reduced 
by 30% (HT' > 200 GeV), 34% (HT' > 300 GeV) and 
39% (HT' > 400 GeV), respectively, compared to no 
cut on HT'. 



We thus can further improve the significances S/ \JB 
and Zq with a cut on HT'. For example, for HT' > 
300 GeV, we obtain S/y/B = 24 and Z a = 8.1. A 
signal is clearly visible. Note, that we cannot give 
meaningful numbers for the significances for HT' > 
400 GeV, because our background samples are not 
large enough. In addition, even if we cut away all 
the backgrounds in Tab. IVIll we still might have some 
small (unknown) QCD backgrounds. Therefore, we 
employ the HT' > 300 GeV cut for our parameter 
scans in the next section in order to get meaningful 
results. It is the main purpose of the HT' > 400 GeV 
cut to show that a (nearly) background free sample is 
possible. 



We conclude, that the cuts presented in Tab. IVIII 
allow a discovery of the benchmark scenario BC1 with 
early LHC data. Scaling down the luminosity, we 
have found that even with an integrated luminosity 
of 200 pb -1 we still get Z > 5! We also want to point 
out that we only used electrons, muons and jets. Our 
analysis does thus not rely on the reconstruction and 
identification of missing energy, &-jets and taus, which 
might be more difficult with early data. In the next 
section we show, that the cuts work also quite well 
beyond BC1. 



B. Discovery Potential with early LHC Data 

We now extend our previous analysis to a more ex- 
tensive parameter region. We still restrict ourselves 
to early LHC data. We focus on regions of the B3 
mSUGRA parameter space with different mass spec- 
tra compared to BC1 but with the same non- vanishing 
lepton number violating coupling, i.e. A121 = 0.032 at 
Mgvt- For that purpose we perform a two dimen- 
sional parameter scan in the -Mi/ 2 _ tan /3 plane around 
BC1 (M 1/2 = 400 GeV, tan/3 = 13). 

We have chosen a scan in M\i 2 and tan/3 for the 
following reasons. Due to the RGE running, every 
sparticle mass increases with increasing M1/2, espe- 
cially the masses of the strongly interacting sparticles 



M M, HI, cf. Fig. p(a)1 and Fig. p(b)1 By varying 



M1/2 we can thus investigate the discovery potential 
as a function of the SUSY mass scale. 

In contrast, changing tan/3 does not affect most of 
the sparticle masses (see Fig. [T5]) and therefore also 
leaves the total sparticle production cross section un- 
changed. 



11(a) where we 



This can be seen in Fig. 
present the total SUSY particle pair production cross 
section (at ^fs = 7TeV) as a function of M 1 / 2 and 
tan/3. For example, increasing M 1 / 2 from 320 GeV to 
500 GeV reduces the cross section from 2 pb to 0.1 pb. 
Note that for M 1/2 = 320 GeV (M 1/2 = 500 GeV) we 
have squark masses mostly around 680 GeV (1 TeV) 
and a gluino mass of roughly 760 GeV (1.1 TeV) ; see 
Fig. |18(b)| and Fig. |18(a)| in App. El But changing 
tan /3 (for fixed Mi/ 2 ) leaves the cross section almost 
unchanged. 

However, increasing tan /3 decreases the stau LSP 
mass, cf. Fig. |18(d)| Increasing tan/3 increases the 
tau Yukawa coupling and thus its (negative) contribu- 
tion to the stau mass from RGE running [26|, [H [5l[ • 
Additionally, a larger value of tan /3 leads normally 
to a larger mixing between the left- and right-handed 
stau j97|. This further reduces the mass of the stau 
LSP. Therefore, with the help of tan (3 we can change 
the kinematics of the LSP decay products and also the 
momentum of the tau from the decay of the Xi into 
the stau (and a tau) . Note that the mass of the Xi is 
nearly independent of tan/3 (see Fig. 18(c) I, because 
it is bino-like around BC1. 

We show the results of the parameter scans in 
Fig. [TTJ The black region corresponds to excluded 
parameter points, with tachyons [3| or which violate 
the Higgs mass or stau mass bounds from LEP [98j |. 
i.e. m h o > 90-114.4 GeV H HJ, dependingon the 
SUSY parameters [H, and m fl > 86 GeV [32]. We 
lowered the Higgs mass bound by 3 GeV to account for 
numerical uncertainties of S0FTSUSY Note > 
that most of the parameter points in Fig. [TT] are also 
consistent with the observed anomalous magnetic mo- 
ment of the muon, with BR(6 — > sj) and with the 
upper bound on BR(B S —> ^ + /i~) [32J. 
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(a)Signal cross section in pb at LHC at y/s = 7 TeV. 




(b) Selection efficiency for the signal events at LHC at 
= 7 TeV. 
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(c) Number of selected signal events assuming an integrated 
luminosity of 1 fb _1 at y/s = 7 TeV. 
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(d) Significance S/ yB assuming an integrated luminosity of 
lfb- 1 at = 7 TeV. 




FIG. 11: Parameter scans in the Mi/2 _ tan(/3) plane. The other mSUGRA parameters are those of BC1, i.e. Mo = Aa — 
GeV and sgn(/i) = +1. Zq is defined in App. iBl 



We present in Fig. 11(b) the selection efficiency for 



the signal events, i.e. the fraction of signal events that 
pass all the cuts in Tab. EID with JJT' > 300 GeV. For 
the scenario BC1 we obtain an efficiency of roughly 
20%. Going beyond BC1, we can see in Fig. 11(b) 



that the fraction of signal events that pass the cuts 
lies mostly between 10% and 23%. We conclude that 
the cuts work also quite well in other regions of the 
stau LSP parameter space than BC1. 
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(a)Estimated contours of minimum required integrated 
luminosity (fb -1 ) to reach S/\/~B > 5 
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(b) Estimated contours of minimum required integrated 
luminosity (fb ) to reach Zrj > 5 assuming a relative 
background uncertainty of 50% 



FIG. 12: Minimum required integrated luminosity for a discovery at y/s = 7 TeV without [Fig. 12(a) and with [Fig. 12(b) 
systematic uncertainties included. The parameters are as in Fig. 1111 



Fig. 11(b) also shows a correlation between the se- 



lection efficiency and tan/3: For fixed M 1 / 2 the effi- 
ciency decreases if tan /3 increases. This behavior can 
be easily understood. As mentioned above, increas- 
ing tan/3 decreases the mass of the stau LSP. In this 
case, the decay products of a light stau LSP have on 
average smaller momenta compared to a heavy one. 
Therefore, less electrons and muons from stau LSP 
decays will pass the cuts in Tab. IVIIl 

This also explains why we very often obtain a bet- 
ter selection efficiency for larger values of M 1 / 2 (and 
fixed tan 0). Increasing M 1 / 2 increases the mass of all 
SUSY particles including the stau LSP. The SM parti- 
cles from cascade decays and LSP decays have then on 
average larger momenta and thus pass the cuts more 
easily. However, when going to very large M 1 / 2 the 
efficiency can again decrease, because the strongly in- 
teracting particles are very heavy and thus their decay 
products are highly boosted. The final state particles 
might then fail the isolation cuts. 

By multiplying the signal cross section [Fig. 11 (a)| 
with the signal efficiency [Fig. ll(b)| , we obtain in 
Fig. 11(c) the number of signal events that pass the 



cuts for an integrated luminosity of 1 fb 1 at <Js = 
7 TeV. For a light spectrum, i.e. M 1/2 < 350 GeV, we 
observe that 100 or more signal events pass the cuts. 
If we go to a heavier spectrum (Mi /2 ~ 500 GeV), 

where the squarks and gluinos masses are ~ 1 TeV, we 
see that the number of signal events reduces to 20 or 
less. Although O(10) signal events are enough to claim 
a discovery of new physics, we do not have enough 
events for a reconstruction of the sparticle masses, es- 
pecially the stau LSP mass. We will address this issue 
in Sec. ESI 

Wc finally present the resulting significances: In 
Fig. 11(d) we give S/\B~ as a naive estimator and 



in Fig. 11(c) Zo as a more realistic estimator. Note 
that Zq is defined in App. [Bl 

With an integrated luminosity of 1 fb -1 at y/s = 



7 TeV, S/yB suggests that stau LSP scenarios up to 
M i/2 ^ 540 GeV can be discovered. Note that M 1/2 = 
540 GeV corresponds to squark masses [gluino masses] 
of 1.1 TeV [1.2 TeV]; see Fig. p(bQl [Fig. p(a)] , in 
App. [A] If we include a 50% systematic uncertainty 
on the background estimate the discovery reach is re- 
duced to M1/2 iS 460 GeV. In this case we have squark 
masses (gluino masses) around 950 GeV (1.1 TeV). We 
want to point out that due to the striking multi-lepton 
signature, the discovery reach for stau LSP scenarios 
with early LHC data is larger than that for i?-parity 
conserving mSUGRA models [85j . 

We have also translated the significances of 



Fig. 11(d) and Fig. 11(e) to the minimum integrated 



luminosity that is required for a five sigma discov- 
ery. For S/\f~B this can easily be done, because S/yB 
scales with the square root of the luminosity. For Zq 
the procedure is more involved, cf. App.[Bj The results 
for S/y/B \Zp\ axe sh own in Fig. p(a)1 [Fig. |12(b)] . We 
can see in Fig. 12(b) that the benchmark scenario BC1 



can be discovered with an integrated luminosity of less 
than 200pb _1 . If the systematic uncertainties are un- 
der good control one might even claim a discovery with 
roughly 50pb _1 , cf. Fig. [12(a) 



C. Stau Mass Estimate 

We have shown in the last section that the stau LSP 
scenario with a non- vanishing A121 coupling can al- 
ready be tested quite stringently with early LHC data. 
If a discovery has been made, one would try to recon- 
struct the sparticle mass spectrum; especially the stau 
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FIG. 13: Invariant mass distribution of the visible part of 
the hardest tau, t , with the two nearest (in AT?) charged 
leptons, I, of the first or second generation. The black line 
(purple line) gives the distribution for the opposite-sign 
(same-sign) lepton pair, £ + £~ (^ ± ^ =F ), plus the tau. The 
distribution is denoted by OS (SS). The green line shows 
the difference of the OS and SS distributions. The red 
histogram corresponds to the correct t^£ + £~ combination, 
i.e. all three leptons stem from the same stau decay. We 
also fitted a Gaussian (blue line) to the green histogram. 

LSP mass. We therefore propose a method to estimate 
the latter. Note, that we do not include systematic 
uncertainties. This has to be done after a discovery. 

In and around BC1, the stau LSP decays via the 
4-body mode 

-tT^t+t-v, (7) 

where £ denotes an electron or muon; see Tab. [TTJ By 
calculating the invariant mass of the stau LSP decay 
products one can in principle reconstruct its mass. But 
due to the initial neutrino and the neutrino in the sub- 
sequent tau decay, not all decay products are visible. 
However, we can still build the invariant mass of the 
£ + £~ pair with the visible part of the (hadronically 
decaying) tau. We then expect a kinematic endpoint 
in the invariant mass distribution, which should lie at 
the true stau LSP mass. 

The black line in Fig. [T3] shows the invariant mass 
distribution of the £ + £~ pair plus the visible part of 
the tau. The distr ibution corresponds to 50 000 BC1 
signal events [l00| . We do not include the SM back- 
grounds, because as shown in the previous section, it 
is possible to select a (nearly) background free data 
sample, cf. Tab. IVIII Leptons are combined as fol- 
lows. In each event we select the hardest identified 
tau. We then look for the two closest opposite sign 
(OS) £s in AR = ^Acj) 2 + Arf. Finally, we calculate 
the invariant mass of these three leptons. If such a 
lepton triplet cannot be found in an event, we discard 
it for mass reconstruction. Note, that we do not em- 
ploy a detector simulation here. We only employ the 
following cuts. Electrons and muons (hadronically de- 
caying taus) arc identified if their (visible) pr is larger 



than 7 GeV (10 GeV). The electrons and muons (tau 
decay products) must also lie within \r/\ < 2.5. We 
use the same setup in what follows if not otherwise 
mentioned. As an additional cut only those selected 
combinations are used, where the distance in AR be- 
tween both leptons and the tau is smaller than 1.5. 
In principle this cut may distort the invariant mass 
spectrum, especially for very high stau masses. We 
checked both options with and without this cut and 
observed, that both give nearly identical results in the 
precision of the estimated stau mass. 

We see in Fig. [151 that the black line has a poorly de- 
fined endpoint at the true stau mass of 148 GeV. This 
is due to combinatorial backgrounds, i.e. we some- 
times combine the wrong tau and electrons or muons 
with each other. Various other combination methods 
have also been studied, for example starting from the 
hardest electron or muon or using A<j> instead of AR. 
A small improvement is possible by vetoing the combi- 
nation r/i + /i _ . This does not increase the fraction of 
correct combinations but reduces the number of wrong 
combinations. The reason is, that the stau LSP can 
not decay to t/jP~ pTv via a coupling A121; cf. Tab. HT1 
We did not to veto the r/i + /i~ combination as this 
keeps the method more model independent. 

In 18% of all events no combination can be found, 
because no hadronic tau in the given kinematic range 
exists. In an independent 18% of all events there is 
only a tau, which does not stem from a stau decay, 
i.e. every method will choose the wrong tau. With our 
method, 30% of the chosen combinations are correct, 
while 25% include a (wrong) tau lepton that stems 
from the Xi decay, which belongs to the same decay 
chain as the chosen leptons. In 14% of the combina- 
tions the tau neither stems from a stau nor a X? decay 
and in 15% at least one lepton stems from the other 
stau decay in the event. In 10% at least one lepton 
comes from another source, i.e. not from a stau decay. 

In order to reduce the combinatorial backgrounds 
and thus to sharpen the kinematic endpoint of the 
invariant mass, we also combine the hardest tau with 
the nearest same sign (SS) lepton pair, £ ± £ ± . The 
respective invariant mass distribution is given by the 
purple line in Fig. [13j We then subtract the r+SS 
distribution (purple line) from the r+OS distribution 
(black line) and obtain the distribution given by the 
green line denoted by OS— SS. 

The OS— SS invariant mass distribution now follows 
much closer the distribution that arises from the cor- 
rect t£ + £~ triplet (red histogram in Fig. IT3"|) . With- 
out the cut on the angular distance between the two 
leptons and the tau the same sign distribution shows 
a long tail at high invariant masses, which leads to 
an over-subtraction at high masses. This can be ex- 
plained by the fact, that one of the same sign leptons 
mostly stems from the other stau decay in the event or 
another source and therefore often has a larger angle 
to the first lepton leading to higher invariant masses. 
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FIG. 14: Stau mass sensitive observable versus true stau mass (see text for definition) for the scenarios presented in 
Sec. IV Bl Different colors and shapes of the points correspond to different intervals of the Xi mass. Fig. 14(a) For 
each parameter point, 10 000 signal events were simulated. Fig. 1 14(b)] Estimates for an integrated luminosity of 5fb _1 
including event selection cuts. Only scenarios where at least 150 events pass the cuts are used. The error bars show 
to what precision the estimated stau mass can be measured. The errors correspond to statistical fluctuations and are 
estimated as described in the text. 



Still one can find stau mass sensitive observables, like 
the intersection of the OS— SS distribution with the 
a;-axis, which show a good correlation with the true 
stau mass. 

We observe that the OS— SS histogram has an end- 
point near the true endpoint. As the stau mass sensi- 
tive observable, we fit a Gaussian on the OS— SS dis- 
tribution and take the value, where it drops to 10% 
of its maximum (marked by a star in Fig. I13[) [101]. 
Although the observable lies below the true stau mass 
of 148 GeV, we can estimate it from the observable as 
long as there is a clear and known correlation between 
the two. This has successfully been demonstrated for 
example in Ref . [EH, [86| . 



This is indeed the case as one can see in Fig. 14(a) 



Here we take the stau LSP scenarios of our parameter 
scans in Sec. IV Bl We then simulated 10 000 signal 
events for each scenario and determined from these 
the estimated stau mass as described above. The dif- 



ferent colors of the points in Fig. 14(a) correspond to 
different Xi masses. 



We can see in Fig. 14(a) a clear correlation between 
the true stau LSP mass and the observable from the 
OS— SS invariant mass distribution. We also see that 
there is only a small systematic dependence of the es- 
timated mass on the \i mass. For example, for a stau 
mass of 120 GeV, the observable can increase roughly 
from 100 GeV to 140 GeV, if the x? mass increases 
from 120 GeV to 240 GeV. This is because a heavier 
X? leads to a harder tau from the x? t\t decay. 



One can use Fig. 14(a) to translate the observable 
to the true stau mass. Such an analysis is even possi- 
ble in a limited way with early LHC data, as can be 
seen in Fig. |14(b)| where we again show the observ- 
able versus the true stau mass. Now, we have only 
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FIG. 15: OS-SS distribution as in Fig.[T3l but here includ- 
ing event selection cuts and randomly selecting events cor- 
responding to an integrated luminosity of 5 fb , i.e. un- 
certainties and fluctuations correspond to those expected 
for 5 fb" 1 . 



included scenarios, where at least 150 events in 5 fb -1 
pass our cuts. Otherwise, we would not have enough 
statistics for the mass reconstruction. We applied the 
event selection cuts that are given in Tab. I VIII with 
HT' > 300 GeV. The error bars correspond to the 
precision with which the observable can be measured 
assuming an integrated luminosity of 5 fb . We do 
not include systematic uncertainties. 

The statistical uncertainties in Fig. |14(b)] were esti- 
mated in the following way. Out of 10 000 simulated 
signal events, events have been randomly chosen to 
get a sub-sample corresponding to an integrated lu- 
minosity of 5 fb -1 . This procedure was repeated 100 
times for each point to obtain different sub-samples. 
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The observable of these sub-samples follows a Gaus- 
sian distribution, where its width corresponds to the 
statistical uncertainties in Fig. |!4(b)| Fig. [TU shows 
an example of one sub-sample for the BC1 scenario. 

We can see that with an integrated luminosity of 
5 fb _1 , rough estimates of the stau LSP mass are pos- 
sible. 



VI. TAU IDENTIFICATION 



We have seen in Sec. IV Cl that the identification (ID) 
of tau leptons is vital for the mass reconstruction of 
the stau LSP. However, as we have pointed out al- 
ready in the discussion of Fig. 2(d) tau ID within the 
framework of BC1 is difficult. We now give a short 
explanation for this observation. 

The ID of hadronic decays of tau leptons in the de- 
tector makes use of special properties of the hadronic 
tau decay. First one expects one or three charged 
pion tracks, i.e. single-prong and three-prong events, 
respectively. Furthermore, jets from tau decays are 
usually rather collimated compared to jets from the 
hadronization of quarks or gluons. The ratio of energy 
deposits in the electromagnetic calorimeter to those in 
the hadronic calorimeter can give further hints. Lcp- 
tonic tau decays (r — > v e v T e and r —> i/^i/ T pi) are usu- 
ally not considered, because the decay length of tau 
leptons is mostly too small to distinguish the decay 
products from prompt electrons or muons. 

However, the special event topologies in BC1 (and 
possibly other bcyond-the-SM models) can make it 
very difficult to reliably identify tau decays in those 
events. In BC1, and also more generally in nearly all 
stau LSP scenarios, a very dense environment is ex- 
pected, i.e. a high multiplicity of charged particles; 
see Tab. [TTJ Therefore, overlaps between the tau jets 
and jets from the SUSY decay chain or leptons from 
the stau LSP decay are very likely. In addition, dif- 
ferent tau jets might also overlap. A reduced tau ID 
efficiency compared to most SM processes is thus ex- 
pected. 

Both Delph.es and PGS [7(| use approaches for 
the detector simulation of tau leptons, which can be 
seen as simplified versions of algorithms for tau lepton 
ID used by the experiments. Observables, like elec- 
tromagnetic collimation and track isolation, are cal- 
culated from the simulated calorimeter deposits and 
tracks. Cuts are then applied to decide, whether a 
jet is tagged as a tau lepton. In both codes we only 
consider 1-prong candidates, i.e. candidates with one 
assigned track. 

We show in Fig. [T7| the efficiency for the ID of 
hadronic tau decays in BC1 and for (SM) Z —> tt + 
ljet production. We also give the fraction of wrongly 
tagged tau jets per event. We show results for the fast 
detector simulation Delphes [Fig. 17(a) Fig. |17(d)| as 
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FIG. 16: Number of hadronically decaying taus in the n- 
Pt plane in arbitrary units. Only the visible parts of the 
taus are considered. No detector simulation and accep- 
tance cuts have been applied. Bins without entries are 
white. 



well as for PGS [Fig. |17(e)f Fig. |17(h)] . In addition, the 
standard overlap removal is applied |l02| . 

Both detector simulations show that the ID effi- 
ciency is significantly lower in BC1 than in Z — > 
tt + ljet. This has to be kept in mind, as many 
methods used in the experiments to estimate the tau 
ID efficiency rely on Z — > tt as a "standard candle" 
|56( . For example, for a (visible) t&u-pr of 45 GeV in 
the inner part of the detector, the tau ID efficiency is 
40% [60%] for Z^tt but only 1 0% [30 %] in B C1 ac- 
cordi ng to D elphes [ PGS], c f. Fig. |17(a)| and Fig. |17(c) 
[Fig. 17(e) and Fig. 17(g)| at the respective working 
points on the cfficicncy-vs-rcjection curve. 

We also see in Fig. [17] that in general the tau ID 
in Delphes has a much lower fake rate than the one 
in PGS. However this comes at the cost of a slightly 
worse identification efficiency giving a different work- 
ing point of the tau ID. As a conservative approach, 
we take the numbers from Delphes. Even in the op- 
timistic case the tau ID efficiency cannot be expected 
to be better than about 25% over a wide pr range 
and can be even lower in the interesting low pr range, 
i.e. (visible) tau-pT < 30 GeV. On can see in Fig. [TBI 
that most tau leptons even have visible transverse mo- 
menta below our ID threshold of 10 GeV. Given that 
only 65% of all tau leptons decay hadronically the ef- 
ficiency with respect to all tau leptons is below 10%. 
This means for less than one in 10 000 BC1 events, 
one can expect all four tau leptons to be correctly 
tagged as such. This explains the low tau multiplicity 



of Fig. 2(d) 



In conclusion the identification of hadronic tau de- 
cays in BC1 and presumably in similar models as well, 
is challenging for the experiments at the LHC. Precise 
estimates of the identification efficiency need studies 
with full-fledged detector simulations. The direct es- 
timate of the efficiency from measured data may be 




(a) ID efficiency in Z 



100 120 140 
MC visible p T [GeV] 

ljet for Delphes 



-S 0.030 
u 

3 0.025 
0.020 
0.015 
0.010 
0.005 
0.000 











A 


lr|l<1.0 












m 


1.0 < Iril < 2.0 
2.0 < Inl < 2.5 












o 








































-A- 














■ 


=)*— i i i i i — i — i 


1 i i i 


■ ■ ■ 


■ 



(b) Fake rate in Z ■ 
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(d) Fake rate in BC 1 for Delphes 
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(c) ID efficiency in Z 
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(f) Fake rate in Z 
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FIG. 17: Tau ID efficiencies and fake rates in BC 1 and Z — > tt + ljet for PGS and Delphes. 
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difficult as the numbers are expected to differ strongly 
between BC1 events and the usual "standard candles" , 
like Z — > tt. 



VII. SUMMARY AND CONCLUSION 

If i?-parity is violated, the lightest supersymmet- 
ric particle (LSP) is unstable. Therefore, it can be 
charged and any supersymmetric particle can be the 
LSP. Within the framework of minimal supergravity 
(mSUGRA), this allows for a large region of the pa- 
rameter space, where the scalar tau (stau) and not the 
lightest neutralino is the LSP. 

We have investigated the LHC phenomenology of 
i?-parity violating mSUGRA with a stau LSP and 
with lcpton number violation. In this model, one 
non-vanishing i?-parity violating operator is present 
at the grand unification scale in addition to the R- 
parity conserving mSUGRA parameters. Focusing on 
pair production of strongly interacting SUSY particles, 
we classified in Sec. IHII all possible LHC signatures at 
parton level; cf. Tabs. lUITVl for an overview. The most 
promising scenarios for the early LHC are those where 
the stau LSP decays via a non- vanishing A121 (or A212) 
coupling. This is because each stau LSP decays purely 
lcptonically via a 4-body decay into two electrons or 
muons, a tau and a neutrino; cf. Tab. HI1 

We have here performed a first comprehensive signal 
over background analysis of stau LSP scenarios at the 
LHC with early data. We employed pair production 
of all SUSY particles as our signal. Our results us- 
ing fast simulations of the ATLAS detector show, that 
the LHC has a good potential to test i?-parity violat- 
ing supersymmetry with a stau LSP and multi-lepton 
final states already in the next two years. The bench- 
mark scenario BC1 (see Tab.|T]) with the 4-body decay 
rjk — > T^i~^l~v has many electrons and muons in its 
final states. This allows the selection of a (nearly) 
background-free sample of O(50) BC1 events using an 
integrated luminosity of 1 fb _1 at -Js = 7 TeV. 

In our event selection used to derive the discov- 
ery potential we avoid using tau-leptons. Naively one 
would expect large numbers of tau leptons in the fi- 
nal states from the decay Xi — > fi + r and the LSP 
decay. However, the reconstruction and identification 
efficiency for the tau leptons is expected to be very 
low, i.e. not larger than 25% (see Sect. IVIl]) . because 
of their small momenta. But a precise estimate is only 
possible with a full detector simulation. Additionally, 
overlaps between jets from tau decay products and 
other jets and leptons in the event further reduce the 
efficiency by a factor 2-3 compared to simple event 
topologies, like Z — > t + t~. Our estimates using fast 
detector simulations predict no identified tau leptons 
in the majority of the BC1 events, even though four 
tau leptons are expected at tree level. Instead, we only 
use the transverse momenta of the leading and sub- 



leading electron and the leading muon, the scalar sum 
of the electron and muon transverse momenta, X^Pt 
and the visible hadronic mass, HT' = X^jct 1-4 Pr- 

In wide ranges of the Mi/2-tan/J parameter 
space around the benchmark point one can achieve 
event selection efficiencies of about 20% for nearly 
background-free samples. Including a systematic un- 
certainty of 50% on the background estimate the dis- 
covery reach is up to M1/2 = 460 GcV for 1 fb" 1 . The 
benchmark scenario BC1 itself can even be discovered 
with an integrated luminosity of less than 200 pb _1 . 

Despite being potentially easy to discover in the first 
LHC data, there arc some difficulties in measuring the 
mass of the stau-LSP in BC1. Due to large combinato- 
rial backgrounds, the selection of correct combinations 
of stau-decay products is difficult and sometimes even 
impossible, because the visible momenta of the tau lep- 
tons are too small to be reconstructed. Neutrinos from 
the stau decay and the successive tau decay only al- 
low for the reconstruction of a smooth endpoint in the 
invariant mass distribution. However, we could show 
that observables can be reconstructed, which show a 
good correlation with the true stau mass, even though 
they will need much more data than 5 fb _1 to provide 
a precise measurement. 
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Appendix A: Properties of the benchmark 
scenario BC1 

We review in this appendix some properties of the 
benchmark scenarios BC1 as described in Ref. [33. 



1. Branching Ratios and Mass Spectrum 

In Tab. IVIIII we show the supersymmetric mass 
spectrum and the branching ratios (BRs) in the BC1 
scenario. The non-vanishing B3 coupling at Mqut is 
Ai2i|gut = 0.032. The other mSUGRA parameters 
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TABLE VIII: SUSY mass spectrum and BRs of the benchmark scenario BC1 [32 
are shown. _R-parity violating decays are in bold face. 



Only decays with a BR of at least 1% 



are M Q = A = GeV, M 1/2 = 400 GeV, tan (3 = 10, 
and sgn(^i) = +1. 

The heavy part of the spectrum, e.g. gluinos and 
squarks, looks very similar to mSUGRA scenarios with 
a Xi LSP, like SPSla [§7j]. In BC1, these sparticles 
mainly decay via two-body decays mediated by the 
usual imparity conserving gauge interactions into a 
lighter sparticlc and a quark. 

Also the middle part of the mass spectrum, i.e. in 
Tab. rVTTTI masses of roughly 300 - 500 GeV, is very 
similar to many Xi LSP mSUGRA scenarios. This in- 
cludes most of the charginos and neutralinos. The 
heaviest neutralinos, xl an( A X4; an d the heaviest 
chargino, xt ' are Higgsino-like whereas the x° an d 
xi are wino-like. These sparticles also mainly decay 
via two-body decays involving gauge interactions. 

However, the light part of the spectrum in Tab. I VILA 
looks very different than the i?-parity conserving 
mSUGRA scenarios with a Xi LSP. We now have a 
stau LSP and the Xi is only the next-to-next-to-next- 
to LSP (NNNLSP). However, it is nearly degenerate in 



mass with the slightly lighter right-handed selectron, 
epi and smuon, fiR. 

We also observe that some of the lighter sparticles 
decay via i?-parity violating interactions, with BRs 
larger than 1%. The stau LSP that can only decay 
via a i?-parity violating operator. Because it does 
not directly couple to the dominant L 1 L 2 £'i-operator, 
it decays in BC1 via a four-body decay, cf. Tab. |TTJ 
The en decays mainly via i?-parity violating inter- 
actions. This is because it couples directly via A121 
and can thus decay via a two-body decay into two 
nearly massless SM particles. The competing i?-parity 
conserving decay mode is a three-body decay, namely 



The Xi NNNLSP is also unstable. 
It decays mainly into a tau plus the stau LSP. See 
Refs. 132, [IH for further details. 



2. Masses and Mass Differences 

We review in Fig. [15] some sparticle masses and mass 
differences in the M 1 / 2 ~t&n (3 plane. The scan is cen- 
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300 350 400 450 500 550 600 300 350 400 450 500 550 600 

M 1/2 [GeV] M 1/2 [GeV] 

(a)Gluino mass. (b)di squark mass. 





300 350 400 450 500 550 600 

Mi/a [GeV] 

(e) Mass difference between the X\ and the fi LSP. 



FIG. 18: Masses and mass differences in the Mi/2-tan/3 plane. The other mSUGRA parameter are that of BC1, 
i.e. Mo = Aq — GeV and sgn(/i) = +1. The blackened out region is excluded due to tachyons or the LEP Higgs or fi 
mass bounds; see also Sec. IV Bl The complete parameter space shown in the figures posses a stau LSP. 



tered around BC1 (M 1/2 = 400 GeV and tan/3 = 13) 
and corresponds to the parameter space for which we 
investigated the discovery potential with early LHC 
data, i.e. Figs. ITT1 and [T2l 



We show in Figs. |18(a)| |18(b)| and |18(c)| the mass 
of the gluino, the left-handed down-type squark, cLl, 
and the Xi, respectively. We observe that these masses 
depend in first approximation only on M 1 / 2 - they in- 
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crease with increasing M 1 / 2 ; see Refs. [5ll. \8i 
ther details. 



for fur- 



However, the stau LSP mass, Fig. |18(d)| is different. 
It also depends strongly on tan /3. For increasing tan (j 
the stau LSP mass decreases. As already described in 
Sec. lVB| this is mainly due to the following effects. On 
the one hand, tan f3 increases the tau Yukawa coupling 
and thus its negative contribution to the RGE running 
of the stau masses. On the other hand, a larger value 
of tan j3 leads in general to a larger mixing between the 
left-handed and right-handed stau. This than results 
in one lighter stau mass eigenstate; cf. e.g. Ref. (27j . 

Therefore, we can change the mass difference be- 
tween the Xi and the stau LSP by changing tan j3 as 



visualized in Fig. 18(e) 



Appendix B: Significance Definitions 

The expected signal significance for a given set of 
event selection cuts depends on the significance def- 
inition used and the estimated systematic uncertain- 
ties. In many phenomenological publications the sim- 
ple S/wB ratio is used, whereas more detailed exper- 
imental analyses use other definitions. We compare 
some definitions in Tab. IIX1 which make different as- 
sumptions on the underlying statistical properties and 
give a very short review of those. For details we refer 
to Refs. [zaU!. 

In all cases the p-valuc of a given number of ob- 
served events under the background-only hypothesis 
with a certain number of expected background events 
is related to the significance Z by 



1 



/2tt 



~ x2 / 2 dx = 1 - 



(Bl) 



where $ is the cumulative distribution of the standard 
Gaussian. 

In all cases we use the so-called Asimov dataset, 
i.e. we estimate the expected number of signal (S) 
and background (B) events by our Monte Carlo esti- 
mate and the expected number of observed events by 
™obs — S + B. Assuming Poisson statistics the p- value 
of observing n bs events when B background events 
are expected is 



oo 



oo 



B" 



n= ^obs n— n bs 
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where Pp is the Poisson distribution. Using Wilks' 
theorem the corresponding significance Zp = — 
pp) can be approximated as 



Z w = \/2(5 + B) ln(l + S/B) - 2S 



(B3) 



in the limit of large statistics. In our case this ap- 
proximation is very good and Zw and Zp have nearly 



identical values. Zyj itself can be approximated with 
the simple S/y/B in the limit of S -C B by expand- 
ing the logarithm. This approximation is clearly not 
valid for our nearly background free sample and is the 
reason for the discrepancy between S/\J~B and Zp in 
Tab. El 

The significance is reduced by uncertainties on the 
background estimate. As the background is estimated 
from Monte Carlo itself, limited statistics of the back- 
ground Monte Carlo translates into uncertainties on 
the estimate of B, which are not included in the pre- 
vious definitions of significance. The ratio between 



Monte Carlo and data luminosity r = 



is used to 



relate the m Monte Carlo background events passing 
the cuts to the background estimate b = m/r. The 
probability to observe n = S + B events is therefore 



P(n,m;S + B,tB) 



Pp{n-S + B) ■ P P (m;rS) 
Pp(n + m;S + B + tB) 
■P B i(n\n + m;p), (B4) 



with p = s +b+ tB arj d Pbi denoting the binomial dis- 
tribution. One can test the background only hypoth- 
esis (5 = 0) with a standard frequentist binomial pa- 
rameter test, giving the p- value 



PBi = 

E p Bi (j\n + m;l/(l + T)) 



(B5) 



The estimated uncertainty of estimating the true br 
by m is \fm. which gives the correspondence Ob = 
\fmjr and with the estimate of b finally r = bja\. 
This relationship can be used to provide an ad-hoc 
background uncertainty o\y to the Binomial p-value, 
cf. Eq. (|B5|) . We used this method in our definition 
of Zq to set a fixed relative background uncertainty of 
/ = 50%, i.e. in summary 



"-off 
n 'on 

Pa 
Z 



l/(b-f) (B6) 

b ■ t' (B7) 

s + b (B8) 

,l on+"off 

E fti(fc|< ff + < n ;l/(l + r'))(B9) 

k—n' on 

s-^i-po) (Bio) 



For large number of background events, as in our sam- 
ple before cuts, this can lead to po > 0.5, i.e. the sig- 
nificance is not well-defined anymore. 

Finally we used a profile likelihood method to incor- 
porate the uncertainty on the Monte Carlo estimate of 
the background in the significance. This has the ad- 
vantage with respect to Z^i , that different scale factors 
Ti for different sub-samples i of the background can be 
treated correctly. A description of the method can be 
found in the documentation of the SigCalc code (89j . 
that was used for the calculation of Zplh- 
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cut 


all SM 
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TABLE IX: Cut flow and different significance definitions for BC1. We assumed an integrated luminosity of J L = lfb -1 at 
yfs = 7 TeV. The given uncertainties include statistical errors only. The significance Zq includes a systematic uncertainty 
of 50%. For the first three rows Zq is not well-defined. 
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because LEP did not explore the possibility of a stau 
LSP decaying via a four-body decay. 

[99] In order to calculate the correct Higgs mass bound 
we used the respective routine implemented in 
SQFTSUSY. 

[100] This corresponds roughly to an integrated luminosity 
of 10 ftT 1 at y/s = 14 TeV. 
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